THE RELATIONSHIP BETWEEN A 
(Q, r) INVENTROY POLICY AND A 
AGE REPLACEMENT MAINTENANCE POLICY 



Thomas Muir Hearn 



United States 
Naval Postgraduate School 




THESIS 



THE RELATIONSHIP BETWEEN A 
(Q , r) INVENTORY POLICY AND AN 
AGE REPLACEMENT MAINTENANCE POLICY 



by 



Thomas Muir Hearn 



April 1970 



This document has been approved public re- 
lease and sale; its distribution is unlimited . 




NAVAL POSTGRADUATE S&HOOIi 
MONTEREY, CALIF. 9394Q 



The Relationship Between a (Q,r) Inventory Policy 

and an 

Age Replacement Maintenance Policy 



by 



Thomas Muir Hearn 

Lieutenant Colonel, United States Marine Corps 
B.S., University of Idaho, 1953 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN OPERATIONS RESEARCH 



from the 

NAVAL POSTGRADUATE SCHOOL 
April 1970 



ABSTRACT 



Cost models are developed to show the relationship 
between inventory and maintenance policies when one component 
of one item of equipment is replaced in accordance with the 
maintenance policy and the components are stocked in accord- 
ance with the inventory policy. A (Q,r) inventory policy and 
failure replacement and age replacement maintenance policies 
are used. The necessary conditions for determining optimum 
values of the reorder quantity, the reorder point, and time 
to replacement are derived. 
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I. 



INTRODUCTION 



The situation to be considered is as follows: There is 

one piece of equipment which will not operate only if a 
certain component has failed. Otherwise it will operate 
satisfactorily. The component fails according to a known 
probability distribution. There is a maintenance policy 
which dictates when the component will be replaced and a 
stock of these components is kept available for use as 
replacements. There is an inventory policy which determines 
the procedures for replenishing the stock of components. 

When the maintenance policy dictates that the component 
is to be replaced, it is replaced if there is a component in 
inventory. If the stock on hand is zero, the equipment be- 
comes inoperative and the time during which the equipment is 
inoperative is called downtime. 

The subject of inventory and inventory policies has been 
examined in considerable depth in many books and papers. In 
this paper only one type of inventory policy will be con- 
sidered, namely a (Q,r) inventory policy. (See Hadley and 
Whitin [1]). This policy is characterized by specifying a 
reorder quantity (Q) and a reorder point (r) . The policy 
requires that an order for a quantity Q of the item be placed 
whenever the stock on hand reaches level r. The length of 
time from the placing of an order to its receipt is called 
the lead time and may be fixed or variable. The demand on 
the inventory system is stochastic and is dependent on the 



probability distribution of failures and the maintenance 
policy . 

A maintenance policy determines when a piece of equip- 
ment will be inspected, repaired or replaced. McCall [2] 
has published a survey of maintenance policies for stoch- 
asticly failing equipment. In what follows, only two types 
of maintenance policies will be considered, failure replace- 
ment and age replacement policies. In a failure replacement 
policy, a component is replaced upon failure; in age re- 
placement the component is replaced upon failure or after it 
has been in use a prescribed length of time, whichever comes 
first. If an item has a continuous, strictly increasing 
failure rate, Barlow and Proschan [3] have shown that there 
exists a unique replacement age which may be infinite, 

(i.e., use failure replacement). 

It is easily seen that if the time to replacement of an 
item is shortened, the quantity used in a given length of 
time will increase. This means that the demand on the 
inventory system will increase. It also demonstrates that 
there is at least one easily seen connection between mainte- 
nance and inventory policies. McCall [2] states that the 
interaction between maintenance and inventory models has 
never been analyzed and suggests that the question merits 
additional research. Falkner [4] considered this inter- 
action in the case in which there is a fixed planning horizon 
and the time to replacement is recalculated after every 
replacement; i.e., a sequential maintenance policy is used. 
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The problem is formulated as a dynamic programming problem 
and methods of determining the optimal values for initial 
inventory and times to replacement are given. 

Although all equipment certainly has a finite life span, 
and optimal sequential replacement policies have been shown 
to be superior to optimal periodic policies for a finite life 
span (see Barlow and Proschan [5]), sequential replacement 
policies often prove to be impractical. Equipment life span 
in the military is usually of indefinite length and is suf- 
ficiently long that it may be considered infinite. The 
application of a sequential maintenance policy requires an 
excessive amount of administrative overhead and requires 
techniques that are usually not known to individuals perform- 
ing maintenance on anything but the most sophisticated 
equipment. For these reasons the infinite horizon models 
considered in this paper are more realistic in many military 
situations. The decision variables are expressed in such a 
way that they can be more easily understood and applied. 

The first case considered is the use of a (Q,r) inventory 
policy and a failure replacement policy. A cost model is 
developed which is used to determine a minimum cost inventory 
policy. For those cases in which cost is not considered to 
be the proper criterion (or sole criterion) for optimization, 
constraints on the length of downtime and the percentage of 
downtime are considered. 

The second case to be considered is the use of a (Q,r) 
inventory policy and an age replacement maintenance policy. 
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A cost model is developed which is similar to the cost model 
for the first case. The model does contain an additional 
decision variable, the length of time to scheduled replace- 
ment of a component. The first case is actually a special 
case of the second case in which the age to replacement is 
infinite. However for clarity we present the two cases 
separately . 
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II. COST MODEL WITH FAILURE REPLACEMENT 



A. DEFINITION OF TERMS 

Reorder Quantity (Q) - The quantity of parts that is 
ordered each time stock is replenished. 

Reorder Point (r) - The inventory level at which the 

order for the reorder quantity is placed. 

Lead Time (T) - The time lag from the time an order is 

placed until it is received. 

Cycle Length (L) - The time from placing of one order to 

the placing of the next. 

Fixed Reorder Cost (A) - The fixed costs of placing an 

order for more inventory. 

Holding Cost (h) - The cost per item per unit time 

charged on all items used during a cycle over the length of 
a cycle. 

Downtime Cost (tt) - The cost of downtime for one unit of 
time . 

Purchase Cost (C) - The cost to buy one item of inventory. 

Failure Replacement Cost (C^) - The cost of replacing 
the failed component on the equipment. This cost is in 
addition to the purchase cost C. 

Length of Downtime (D) - The length of time in one cycle 
during which the equipment cannot be operated. 
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B. SITUATION AND ASSUMPTIONS 



Consider the situation as described in the introduction, 
in which a piece of equipment is maintained under a failure 
replacement policy and a (Q,r) inventory policy is used to 
supply the replacement parts. The lifetimes of the components 
are independent, identically distributed random variables 
with a known distribution function, F(x) . This function is 
assumed to have an increasing failure rate, that is 
f(t) . . 

— - — — is increasing in t, 

F(t) 

where f(t) is the density function and F(t)=l-F(t). The 
probability of failure of an item with an increasing failure 
rate increases with time in service, i.e., it exhibits wear- 
out . 

In the development of the model further assumptions are 
made. It is assumed that the time required for component 
replacement is zero, thus downtime occurs only when the stock 
of replacement parts is exhausted. The lead time is assumed 
to be fixed. The cost of downtime is assumed to be constant 
throughout the length of the downtime. Figure 1 illustrates 
the inventory level in a system using a (Q,r) inventory 
policy and a failure replacement maintenance policy. It is 
assumed that Q>_ r+1, this insures that when an order is 
received, the quantity on hand is at least r. 
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Inventory Levels with Failure Replacement Policy 



C. MODEL FORMULATION 



In formulating the model, an expression is first devel- 
oped for the costs for one cycle; these costs are then 
averaged over an infinite horizon. 

1 . Costs for One Cycle 

The costs for one cycle are expressed as; 

cost of placing one order + purchase cost + holding 
costs + costs of downtime + cost of replacing the component. 

Thus, the total cost, K , for one cycle is: 

K=A+CQ+hQL+TD+C f Q (1) 

In the above equation, the quantities L, cycle 
length, and D, length of downtime, are random variables. 
Expressions for the expectations of these random variables 
are now found. 

2 . Expected Length of Downtime 

As previously stated, downtime will occur if the 
component in service fails when the inventory level is zero. 
Since the inventory level is r when an order is placed, r+1 
failures must occur prior to the end of the leadtime T for 
downtime to occur. An example of this can be seen in 
figure 1. 

Let = time to failure of the ith component, 

Y = length of downtime in a cycle. 

Since downtime cannot be less than zero and the 
leadtime T is a fixed quantity, then 
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r+1 

Y = Max [0 ,T- Z X i ] . 



i=l 



n 

Defining S = Z X. 

n i=l 1 



then Y = Max [0 ,T-S r+ ^ ] . 



Let the distribution of Y be G(y) = P [Y<y] . Then 
G (y ) = P[S r+1 >T-y]=F r+1 (T-y) y<T , 



( 2 ) 



where F (x) is the n-th convolution of F(x) with 



n 



itself. 

Let 6 = E [ Y ] = expected length of downtime per cycle, 

6 = E[Y] = /J [l-G(y) ] dy = /J [1-F r+1 (T-y) ] dy = /J F r+1 (x)dx. 

3. Expected Length of a Cycle . ^ 

A cycle has been defined as the length of time from 
the placing of one order to the placing of the next order. 
During a cycle there are exactly Q failures, each preceded 
by a length of time X^ , during which the equipment is 
operating . 

During a cycle the length of time the equipment is 

Q 

operating is Z X. , the length of time the equipment is not 
i=l 1 

operating is D, the downtime. Thus the length of a cycle is 



Q 

L = Z X.+D. 



i-1 1 



Taking expectations and letting £. = E[L] 



Q 

1 = E [L] = Z E [X. ]+E [D] . 
i=l 1 
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Let E[X]=y the mean time between failures (MTBF) . 



Thus, 1 = Qy + /q F r+1 (x)dx . 

4 . Costs for an Infinite Horizon 

From equation (1) , the average costs per cycle may 
be determined by taking expected values. That is, the costs 
are averaged over an infinite horizon. Letting K = the cost 
for one cycle, 

E [K] = A + CQ + hQE [L] + ttE [D] + C f Q 
= A + CQ + hQ£+ it 6 + C f Q . 

The sequence of cycles is a renewal process since the 
cycle lengths are independent, identically distributed random 
variables. Let the expected number of renewals in [0,t] be 
H(t). The expected cost for the interval [0,t] will then be 
E [Cost] = H (t) E [K] + e (t) , where e ( t) /t 0 as t-* 00 , 
since e(t) represents the cost in the part of a cycle in 
[0,t] and hence is bounded. 

The average cost per unit time is 



lim 

t+0 



E [Cost] 
t 



lim 

t+0 



E [K] + ] 



=E [K] 



limH ( t ) 
t+0 t 



+ 0 . 



time is : 



• H ( t ) 1 

Since lim — = — 

t-0 1 * 



the average cost per unit 



% = 



E [K] 

1 



A + CQ + C f Q + tt 6 + hQ 

l 



(4) 
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Since X is a function of Q and r equation (4) may 



be written 

A +(C+C f )Q +T r /JF F (x)dx 

K(Q, r) = ±— - ^ +hQ . 

Qy+/ 0 F r+1 (x)dx 

5 . Properties of Cost Function 



( 5 ) 



T 

6(r) = /q F r+ -^(x)dx is a decreasing function of r. 

It follows from this that ^{(Q,r) is either monotone increas- 
ing or monotone decreasing in r, depending on the relative 
magnitudes of certain costs in equation (5) . Let C'=C^+C. 



A+C 1 0 

If ^ ' Qy > b K is monotone increasing in r. 



A+C 1 Q 

If — rcQ - ” — < b K + s monotone decreasing in r. 

Since -1<_ r £Q-1 the optimum value of r is 
, . A+C 'Q . , 

r = - 1 lf ~ToiT > 1 ' 

r = Q-l if < 1 . 

ttQm 



( 6 ) 



A+C 1 Q 

If ■ = 1 the function is constant in r. 

At r = -1, 6=T, thus 



% (Q/ -1) 



A+C 'Q+ttT 
Qy+T 



+ hQ . 



Assuming K is continuous and differentiable in 
Q the optimum value of Q is the solution to d'K/dQ=0 . From 



the above , 



dX C’ A+C ' Q+ttT , 

dQ - Qp+T - (Qy+T) 2 

after solving for Q , 
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_ ) Ay + (fry -C 1 ) T T 

y v 2 y 

hU 



For large values of r, 6-0, thus 



K(Q,r) = + hQ * 



In the same manner as above, the optimum Q for large 
values of r is: 



Q - / |r • 



(7) 



This is nearly identical to the well known lot- 
formula [1] and would be identical if the assumption had 
been made that holding costs were proportional to Q/2 . 

For large values of r, the second derivative of K 
with respect to Q is 



d K _ 2y (yA+yC’Q) _ 2yC' 



dQ z 



(Qy) 



(Qy) 



2A 

Q 3 y 



> 0 



Thus X(Q, r) is convex in Q for large r. 
the second derivative of 



is 



dQ 2 



_ 2y (yA+yc 'Q+ u yT) _ 2yc' 



For r =-l , 



(Qy+T) J 

2y 2 A+2 -fru 2 T-2yC’T 



(Qy+T) 



(Qy+T) z 

This will be positive if y (A+irT) > C'T which is then 
the necessary condition to have convexity in Q. 
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D. DOWNTIME RESTRICTIONS 



It is frequently desirable, particularly in the military, 
to maintain the length of downtime or the percentage of 
downtime below a certain level. The necessary conditions for 
applying these restrictions have been developed and the 
problem of applying them subject to cost constraints has been . 
considered . 

1 . Restrictions on the Length of Downtime 

If D = downtime then, from equation (2), 

P [ D <y ] = F r+1 [T-y] 0<y<T 

= 1 y >T 

Note that the length of downtime is a function of r 
alone, given T as fixed. If it is desired to control the 
length of downtime this can be done by varying r. For 
example, since downtime will always be a fraction of T, to 
keep this fraction less than or equal to a, with probability 
at least 3, we must choose a and 3 such that P[Y<aT]>_3. 

That is, 1-3 > F r+ i [(l-a)T]. 

Since P r+ ^ (x) is a decreasing function of r and 

■ L;Lm F , , (x) = 0, there is a minimum value of r for which the 
r->-°° r+ 1 

above inequality holds. 

2 . Fraction of Downtime 

The expected fraction of time during which the item 
of equipment will be down is: 
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expected length of downtime 

fraction downtime = 

expected length of cycle 



6 

Qy+<5 



/ o F r+i ( x ) dx 

= . We wish to investigate 

Qy+/ o F r+1 ( x ) dx 



how this function varies with r and Q. 

The following lemma will be needed: If A>0 and B>0, 

then BTAigfSr £or e>0 - 

Proof: assume 



A < - A - C , then A ( A+B) -As <A ( A+B) - Ae-Bc and 



B+A B+A-e 

0 < -Be, but this contradicts the fact that e>0 and B>0, 
therefore, ^ A , > A e 



B+A - B+A-e ' 

The maximum amount of downtime that can occur is T, 
This will occur when r= -1. For any r>-l, 6<T. 

Thus, from the above lemma 

T 



fraction downtime < 



Qy+T 



It can be seen from the above that by increasing Q 
the fraction of downtime may be made as small as desired. 
Note also that 6 is a decreasing function of r, so, subject 
to the constraint Q>r, by increasing r, the fraction down- 
time may be made as small as desired. 

3 . Cost Constraints on Downtime 

By increasing Q and r and keeping Q>r both the 
expected length of downtime and the fraction of downtime may 
be made as small as desired. Normally there will be a cost 
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constraint of some type; the problem can then be expressed 



as 



Minimize (fraction downtime) = 

subject to K(Q/ r )£k 

where b is a cost constraint; 



6 

QU+6 



or. 

Minimize K(Q/^) 

subject to Qjj+ 5 * <P* 0<p<l 

T 

in this case if <p the constraint is not active and 

need not be considered in the minimization. 

If the length of downtime is of concern the problem 

is. 

Minimize X(Q/ r ) 

subject to 1-g >_ F r+1 [(l-a)T] 
or 

Minimize F r+ ^[(l-a)T] 
subject to ^(Q,r)<b . 
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III. COST MODEL WITH AGE REPLACEMENT 



The cost model with age replacement is a generalization 
of the model with failure replacement presented in II. 

A. DEFINITION OF TERMS 

The following definitions are required in addition to 
those given in II. A. 

Time to replacement (t ) - The time between the placing 
of a component in service and its scheduled removal assuming 
the component does not fail before t . 

Scheduled Replacement Cost (C s ) - The cost of replacing 
a component at the scheduled time to replacement (t ) . 

B. SITUATION AND ASSUMPTIONS 

The situation is the same as previously described. The 
same assumptions still hold. In addition, it is assumed that 
the cost of scheduled replacement of a component is less than 
the cost of failure replacement, C g <C^. 

C. MODEL FORMULATION 

The development of a cost model for an age replacement 
policy follows very closely the development in II. There are 
some additional costs to consider and in addition to the 
distribution of failures F(x), the distribution of time 
between replacements must be known. 
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1 . 



Distribution of Time Between Replacements 



The use of an age replacement policy means that an 

item in service will be replaced upon failure or if the time 

in service reaches t , whichever occurs first. The distri- 

o 

bution of time between failures is still F(x), a continuous 

A 

function with an increasing failure rate. Defining F(x) = 
distribution of time between replacements with an age 
replacement policy, F(x) has the following distribution: 

F (x) = F(x) °£ x<t 0 (8) 

F(x) = 1 . x>t 

— o 

Let y be mean time between replacements, then 
- t o 

U = /° [l-F(x) ] dx . 

0 

2 . Expected Costs for One Cycle 

The expected costs for one cycle are: 

ordering cost + purchase cost + expected holding 

costs + expected costs of downtime + expected costs of 

replacing items upon failure + expected costs of replacing 

items upon schedule at time t . 

o 

If X is the time to failure of a component, then 
the probability of the item failing before the scheduled 
replacement is P [X<t ] = F(t Q ) . For one cycle, the expected 
cost of replacing items upon failure is the cost of a failure 
replacement times the expected number of items that will 
fail, which is the number of items used in a cycle, Q, times 
the probability of failure, F(t Q ). Thus, the expected cost 
of replacing items upon failure = C^QF(t ). Similarly, 

the expected cost of replacing items upon schedule is 
C g Q[l-F(t 0 ) } . 
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The expected costs for one cycle. = 



A + CQ + hQt + tt< 5 + c f QF(t o ) + C s Q(l-F(t o )) . 

3 . Costs for an Infinite Horizon 

Using the same argument as in II. C. 4 the costs per 
unit time averaged over an infinite horizon are obtained. 

A + CQ + C,QF(t ) + C Q ( 1-F ( t ) ) + TT 6 

% = i- 2_ 5 2 +hQ . 

t 

A A 

Recalling that £ and are functions of t , this 

may be written 



ft (Q r F f t Q ) 



T ^ 

A+Q(C+C f )F(t o )+Q(C+C s )F(t o )+7r/ 0 F r+1 (x) dx 



Qy + f Q F r+1 (x)dx 



+hQ . 



4 . Properties of the Cost Function 

a. rjn /s 

Since 6(r) = /q F^ + ^(x)dx is a decreasing function 
of r, the same type of argument as in II. C. 5 gives the 
following for optimum values of r; 

letting C'=C+C^ as before and letting C''=C+C s 
A+C ' QF ( t ) +C' 'QF(t ) 

r=-l if — >1 . 

iTQy 

A+C 1 QF ( t ) +C ' ' QF ( t ) 

r=Q-l if <i . 

iTQy 

Note that the optimum value of r is dependent on 
both Q and t Q . 

/s. 

At r=-l 6 = T, using this fact and assuming X is 
continuous in Q, by taking partial derivatives of 
X(Q,-l,t o ) with respect to R and t Q the required conditions 
for optimum Q and t Q can be found. 
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C ' F ( t )+C''F(t ) A+QC ' F ( t ) +QC ' ' F ( t ) +ttT 
o A _ v, i ° ° o o 

3Q ,„ A , mV 2 



Qy+T 



(Qy+T) 



setting this equal to zero, we find 



Q = 



y (A+ttT) - ( C ' F ( t ) +C ' ' F ( t Q ) ) T 



hy‘ 



T 

y 



( 9 ) 



Taking partial derivatives with respect to t 



„ 4 , QC'f(t ) -QC ' ' f ( t ) ~ i A+QC ' F ( t ) +QC ' ' F ( t )+ T 

= 2 2_ _ Ou 2 ° 

3t ~ ~ 2 

o Qy+T (Qy+T) 



note that M=/q F(x)dx 



and 



y' = 3y/8t Q = F(t ) by Leibnitz rule 



Q ( C ' -C ' 1 ) f ( t ) A+QC ' F ( t )+QC''F(t )+ T 

° = * “ - QF(t o } 3—2 2 



Qy+T 



(Qy+T) 



The optimal t is the solution to the above equation. 

Providing the necessary conditions for r=-l hold 
throughout the region under consideration, an iterative 
procedure may be used to determine the optimum Q and t Q . 

Assuming that r =Q-1 and Q is sufficiently large, 
the necessary conditions for a minimum ft(Q,r,t Q ) are: 



Q = 



hy 



and 



Q (C ' -C ' ' ) f ( t ) F ( t ) 



Qy 



Qy 



2^- (A+QC'F(t o )+QC' ’F(t o ) )= 0 . 
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Rearranging this equation we get 



~ ( to) c ' 1 A 

M F(t } ' F(t o } = C'-C' ' + Q (C ' -C ' ') • 



This is very similar to the equation given in refer- 
ence [2] for determining the optimum t Q when an age replace- 
ment policy is used. The above equation differs only by the 

A 

inclusion of the term ) " on the right hand side. If 

A is zero, it is identical to the equation in [2] which does 
not consider costs due to purchasing and holding inventory. 
Note also that for large Q or large C the effect of this term 
is small. 

Following the same procedures as in II. C. 5 it can be 
shown that ft(Q,r,t Q ) is convex for large values of r and that 
for r=-l the following condition must hold for convexity in 
Q, 



y (A+tt) > (C'F(t ) +C' ' F ( t ) ) T 



D. DOWNTIME RESTRICTIONS 
1 . Length of Downtime 

Using the same reasoning as in II.D.l but with the 
distribution of time to replacement given by (8), to keep 
the downtime less than or equal to a fraction a of T, the 
lead time, with probability 8; the following must hold. 

H> F r+1 [(l-a)T] . 

For any r, and letting k= 1-a , 
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F r+l (kT) = 1 * kT > ( r+1)t 0 

= F r+1 (kT) . kT <(r+l) t Q 

Thus, if (l-a)T > (r+l)t o , 3 must equal zero, that 
is if aT < T-(r+l)t then 3=0. This means that the lower 
limit on downtime is T- (r+l)t . 

setting T-(r+l)t less than or equal to zero 

we t o > 

T 

This implies that if t < — j-=- downtime will occur 
every cycle and it will be no less than T- (r+l)t Q . 

This result indicates that a strict age replacement 
policy is defective in some manner. This is the case and 
this defect will be discussed further in V, where the 
effects of a more realistic policy will be examined. 

For ( 1-a) T< (r+1) t , 

1 - 3 *F , , [(l-a)T] must hold and as in 
— r+1 

II.D.l there is an r' such that 
for all r>r' 

l-3>F r+1 [ (l-ct)T] for all a and 3 • 

2 . Fraction of Downtime 

The expected fraction of downtime is 



5 



Qy+<5 



/ F 
J 0 r+1 



/\ T* 



(x) dx 
F r+1 ( x) dx 
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in a manner similar to II. D. 2 it can be shown that 



. T 

fraction downtime^ — — . 

Qy+T 

/\ /N 

Since 6 is nonincreasing in t and y is increasing 

in t Q , the fraction of downtime is decreasing in t . This 

follows from the above lemma. The fraction of downtime will 

reach a minimum with respect to t at t = °°, i.e., when a 

o o 

failure replacement policy is used. 

This means that the fraction of downtime can be 
decreased by increasing Q or r as was seen in II. D. 2 and in 
addition a decrease in downtime can be brought about by 
increasing t . 

3. Cost Constraints on Downtime 

Constrained minimization problems similar to those 
given in II. D. 3 can easily be developed if required. 
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IV. NUMERICAL EXAMPLES 



A. INCREASING VALUES OF r 

Since the optimal solutions for large values of r 
obtained in II. C. 5 and III.C.4 are based on the assumption 
that the expected value of downtime is close to zero, it is 
desirable to see how rapidly this value approaches zero as 
r increases. Figure 2 contains graphs of the expected value 
of downtime for increasing r for two different failure 
distributions with the same mean. The values are graphed 
for different values of T. Note that in all cases for r 
greater than 7 the value is negligible. 

B. DETERMINING OPTIMA 

The following values are used in this example: 
A=$5.00,TT=$50.00,C=$2.00,C f =$1.00,C s =$0.25, 
h=$0.10,T=3 months or .25 year,y=l month or .0833 years. 

The failure distribution of the component being replaced 
is uniform. 

The optimum value of r will be Q-l for all Q greater 
than or equal to 5 if a failure replacement policy is used. 
This follows from equation (6) . The optimum value of Q was 
determined from (7) to be 24.40. A computer routine that 
searched for the Q that gave a minimum cost found an optimum 
Q of 25 and a cost of $42.90 for a failure replacement 
policy. 
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Exptcttd Downtime Expected Downtime 



.50 



Uniform Distribution 
H = 0833 




r 




r 

Figure 2 
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Using the computer search routine the optimum values of 
Q,r, and t were determined to be Q=25 , r=24, t o =0.145 years 
or about 53 days with a cost of $42.26. 

C. SENSITIVITY ANALYSIS 

Using the values given above a sensitivity analysis was 
conducted to determine the effect of changes in the para- 
meters and the variables. Using the optimal Q of 25 the 
following costs are obtained for the values of t shown: 



t 

o 


Cost 


.08 


50.039 


.09 


47.736 


.10 


45.248 


.11 


44.084 


.12 


43.238 


.13 


42.474 


.14 


42.260 


.15 


42.343 


.16 


42.659 


oo 


42.896 


* 





As can be seen from the table it is better to have a 
failure replacement policy than to replace components too 
frequently. A failure replacement policy entails a cost 
increase of 1.2 percent whereas monthly replacement 
(t Q =.0833) will increase the cost by approximately 18 percent. 

The cost increase caused by varying Q plus or minus 5 
while t Q was held constant at .145 was less than 0.5 percent, 
indicating relative insensitivity to changes in Q. 

The following table shows costs and downtime percentages 
at the optimum values of Q and t when r is varied through 
the range of values for which it has an effect on the model. 
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r 


Cost 


Percent 

Downtime 


-i 


43.43 


11.3 


0 


43.07 


7.8 


1 


42.69 


4.2 


2 


42.43 


1.6 


3 


42 . 31 


0.4 


4 


42.28 


0.1 


5 


42.27 


0.02 


6 


42.27 


0.00 


7 


42.26 


0.00 



The effect of varying the parameters in the model was as 
follows : 

Varying h, the holding cost, had a significant effect on 
the optimum Q and a small effect on the optimum t . The 
changes in cost were significant. 

Varying A, the fixed reorder cost, affected the optimum 
Q and the cost but did not have any sigificant effect on the 
optimal t Q . 

The cost K varies inversely with changes in y, the mean 
time between failures and there are significant changes in 
optimal t . The optimal Q varies by only a small amount as 
y changes. 

Varying tt , the cost of downtime, has no effect until tt 
is low enough that the optimal value of r is -1. 

Varying the costs, C,C^ and C g has a direct effect on ft 
and causes a significant change in the optimum t . The 
optimum Q changes only slightly. 
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V. AN IMPROVED AGE REPLACEMENT POLICY 



A- DESCRIPTION OF POLICY 

The age replacement policy under which the model in III 
was developed requires that a component be replaced upon 
failure or when the time in service reaches t . This policy 
is followed regardless of the inventory level, even if the 
inventory level is zero. This means that the component in 
service will be removed and the equipment will go down. As 
was seen in III.D.l this can lead to a situation in which 
downtime will occur every cycle. 

It is not reasonable to expect that a serviceable 
component will be removed forcing a piece of equipment to go 
down when there is no replacement available. For this reason 
the following improved policy is set forth: 

Follow the age replacement policy as set forth in III 

unless the inventory level is zero. In that case remove the 

component in service only upon failure or upon receipt of 

more inventory if the component has more time in service 

than t . 
o 

B. MODEL FORMULATION 

Figure 3 depicts the inventory levels when such a policy 
is followed. Comparing this policy with the previous policy, 
the following can be noted: 

(1) Cycle lengths will be identical if r is greater 

than -1. 
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Invtntory Levels with Improved Age Replacement Policy 
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(2) The distribution of downtime will be different. 

(3) The distribution of time to replacement of one 
of the components used in a cycle will be different if the 
inventory level reaches zero. 

1. Identity of Cycle Lengths 

If the inventory level does not reach zero or if the 
level reaches zero, but the age of the item in service does 
not reach t , the length of a cycle will be the same under 
either policy. 

If the inventory level is zero and the age of the 
item in service reaches t before more inventory is received 
the two policies require different actions: under the 

previous policy the item in service will be removed, downtime 
will begin and last until a new component is received; under 
the improved policy the item will be left in service until 
failure or receipt of a new component. Note that in both 
cases the cycle ends when more inventory is received; thus 
the cycle lengths will be identical although the lengths 
of downtime will differ. 

Since the cycle lengths are identical to the case 
previously considered, the distributions of cycle lengths 
are the same and the expected length of a cycle (&) will be 
the same. 

If r=-l and the item in service when the inventory 
level is zero is not removed until it fails, the cycle length 
will be greater than or equal to the cycle length in the 
previous case. This case will be considered separately. 
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2 . 



Distribution of Downtime 

Let Y = life of item in service when inventory level 
is zero. 

Y is distributed as F(x) if downtime occurs. 

Letting D = downtime for a cycle 



D = Max [ 0 ,T-S -Y] , 

P [D<t] = P [T-S -Y<t ] 

= P [S +Y>T-t ] . 

Letting G(t)= F *F(t), where the * denotes 
convolution 

P [D<t] = G(T-t) , 0<t<T 

P [D> t ] = G(T-t) . 0<t<T 



The expected length of downtime using this policy is 

then : 

~ T T 

6 ' =E [D] = / G (T-t) dt = / G (x) dx . 

Note that lim 6'=0 . 

r- >°° 

3 . Expected Length of a Cycle 

As a result of V.B.l the expected length of a cycle 

/\ /\ 

must equal Qy + 6 . 

Using the improved policy an expected cycle length 
using the same reasoning as in II. C. 2 is 

(Q-l) y + E [Y] + 6 ' . 
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Equating the two cycle lengths; 



A A 



Qy + 6 = (Q-l)y + E[Y] + 6' , 



AAA 



E [Y] = y + 6 - 6 ' 



= y + / Q F r+1 (t)dt - J Q F r *F(t)dt 



= y + f Q F r *(F(t)-F(t))dt . 



Note that; lim E[Y] = y . 
r — >°° 

The expected length of a cycle can be expressed as: 



Qy + f Q F r *(F(t)-F(t))dt + 6' 



or Qy + 6 . 



4 . Cost Model 

In the same manner as in III.C.3 we can develop the 
following cost model;. 



/\ 

A+C 1 QF ( t )+C''QF(t ) +7T 6 1 

X(Q,r,t ) = ■ 2 +hQ . 

Qy+6 



5 . Properties of the Cost Model 

As in the previous models ft is either monotone 
increasing or monotone decreasing in r, so the optimum r 
will be either -1 or Q-l. 

For large r the model is identical to the model in 
III.D.4 and the necessary conditions for optimum Q and t Q 
are the same. 
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For r=-l, 6' = T and the expected length of a cycle 

A 

will be: (Q-l) y + y + T . 

As mentioned previously, this is not identical to 
the cycle length under the previous policy, a situation that 
can arise only if r=-l. 

Under this expected length of a cycle and the 
expected length of downtime from 2, we can determine an 
optimum value of Q by assuming continuity in Q. It is 



Q = 



y (A+ttT) - (C'F(t ) +C' 'F(t ) ) .(T+y-y) 

— 

hy z 



y-y+T 

A 

y 



Since y-y is positive, Q will be less than in 
equation ( 9) . 

The optimal t will be the solution to the following 
equation : 



0 



Q (C ' -C ' ' ) f (t ) A+QC ' F ( t )+QC"F(t )+itT 

— 7 7T- 2 - - QF(t ) r- 2 t ~2 

Qy+T+ (y-y) (Qy+T+(y-y)) 



An iterative solution of these two equations will 
yield the optimal Q and t for the improved policy when the 
optimum r is -1. 
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VI. SUMMARY AND CONCLUSIONS 



A. SUMMARY 

Cost models are developed to show the relationship 
between inventory and maintenance policies when one component 
of one item of equipment is replaced in accordance with the 
maintenance policy and the components are stocked according 
to the inventory policy. The necessary conditions for 
determining optimum values of Q, the reorder quantity; r, 
the reorder point; and t , the time to replacement are 
derived. 

The models are monotonic increasing or decreasing in r, 

thus the optimal r is either -1 or Q-l. If Q is large 

(approximately 8) further changes in r have little effect, 

thus r could be considered constant. When r is large the 

necessary conditions for optimal Q and t are very similar to 

well-known conditions for determining optimal values of Q 

and t . 
o 

B. CONCLUSIONS 

The cost models can be used to determine optimal values 
of Q, r, and t Q . The optimal values of Q and t are not very 
sensitive to changes in each other or changes in the para- 
meters. The optimum value of r, having only two values, is 
very sensitive to changes in Q or t Q or the parameters when 
near the boundary between increasing and decreasing f(. The 
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cost, however, does not appear to vary much even with changes 



in r and it might be best to use an r large enough 
strain downtime in a specified manner when in this 



to con- 
region . 
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